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ABSTRACT 



Context. The Ultra-luminous infrared galaxy Mrk 231, showing signs of both black hole accretion and star formation, exhibits very 
strong water rotational lines between A = 200 - 670 lum, comparable to the strength of the CO rotational lines. High redshift quasars 
also show similar CO and H2O line properties, while starburst galaxies, such as M82, lack these very strong H2O lines in the same 
wavelength range, but do show strong CO lines. 

Aims. We explore the possibility of enhancing the gas phase HtO abundance in X-ray exposed environments, using bare interstellar 
carbonaceous dust grains as a catalyst. Cloud-cloud collisions cause C and J shocks, and strip the grains of their ice layers. The 
internal UV field created by X-rays from the accreting black hole does not allow to reform the ice. 

Methods. We determine formation rates of both OH and H2O on dust grains, having temperature Tjust = 10 - 60 K, using both Monte 
Carlo as well as rate equation method simulations. The acquired formation rates are added to our X-ray chemistry code, that allows 
us to calculate the thermal and chemical structure of the interstellar medium near an active galactic nucleus. 

Results. We derive analytic expressions for the formation of OH and H2O on bare dust grains as a catalyst. Oxygen atoms arriving 
on the dust are released into the gas phase under the form of OH and H2O. The efficiencies of this conversion due to the chemistry 
occurring on dust are of order 30 percent for oxygen converted into OH and 60 percent for oxygen converted into H2O between 
Tdust = 15-40 K. At higher temperatures, the efficiencies rapidly decline. When the gas is mostly atomic, molecule formation on 
dust is dominant over the gas-phase route, which is then quenched by the low H2 abundance. Here, it is possible to enhance the warm 
(T > 200 K) water abundance by an order of magnitude in X-ray exposed environments. This helps to explain the observed bright 
water lines in nearby and high-redshift ULIRGs and Quasars. 

Key words. Galaxies: starburst, AGN; ISM: chemistry 
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1. Introduction 

. Warm molecular gas (T > 100 K) cools mainly through CO, 
' H2O, and H2 ro- vibratio nal transitions jNeufeld & KaufrnanI 
' ll993HNeufe1d et al.llT995l) . The ground state lines of CO, up to 
• J - 6-5, are extensively s tudied from the ground in a large 
' number of galaxies (see, e.g.. lPapadopoulos et al. l l2OT0h up to a 
redshift z ~ 0.04. Unfortunately, the atmosphere makes it im- 
possible to study CO transitions with 7 > 8 and H2O rotational 
lines. It is these lines, however, that are tracing the warm molec- 
ular gas. Moreover, these lines are sensitive to extreme UV ir- 
radiation from star-formation and X-rays due to supermassive 
black-hole accretion. Some high-7 CO lines have been observed 
in bright galactic objects such as the Orion Bar and Sagittarius 
A (cf. Hollenbach & Tielens 1999), and the Class object 
L1448-mm (Nisini et al.J999.) . It is only recently that significant 
progress is being made in the study of the warm molecular in- 
terstellar medium (ISM) in nearby galaxies. Both high-7 as well 
as rotational water lines have been observed abundantly after the 
launch of the Herschel Space observatory, that is operating be- 
tween wavelengths A = 55-672 fim. The successful HerCULES 
open time key program has observed 29 close (U)LIRGs with 
both SPIRE and PACS instruments, revealing CO lines up CO 
J = 14-13 (while current follow-up PACS observations 
might detect even higher transition) and H2O lines in Mrk 
23 1 that are of comparable strengths dvan derWerf etal1l2010t 



iGonzalez-Alfonso et aI.ll20Toh . These CO and H2O lines are at 
least in part produced in regions of the galaxy exposed to X- 
rays, so-called XDRs (Maloney et al. 1996; Meijerink & SpaansI 
2005), but emission produced through shock e xcitation may 
also contribute as observed for OH in Mrk 231 (iFischer et al.l 
2010). Moreover, Herschel-PACS also detected a lot of high-7 
CO lines as well as H2O and OH in e mbedded objects, such as 
DK Chamaeleontis dvan Kempen et aL,2010) and in disk atmo- 
spheres, e.g., HD 100546 dSturm et al.ll20'Tol) . where X-rays also 
affect the chemical-thermal structure. Even though this paper fo- 
cuses on an application to ULIRGs, it will have applications to 
other objects as well. 

Understanding the excitation and chemistry of these impor- 
tant molecules is therefore key in studying the dominant phys- 
ical processes in the ISM and star-formation. CO is a linear 
molecule, and the radiation transfer models allow to calculate 
far-infrared to millimeter spectra relatively easy for a given 
physical structure. Models predicting CO emission, where both 
the thermal-chemical structure and radiation transfer are calcu- 
lated simultaneously have more problems to reach agreement. 
The case for H2O is even worse, due to its complex nature: high 
critical densities ncrit > 10^ - 10'" cm""* in combination with 
high opacities, infrared-pumping and maser actively, makes an 
excitation calculation extremely challenging. Also, the excita- 
tion state in which water will enter the gas-phase after desorp- 
tion or evaporation from a dust grain is very uncertain. A first 
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guess in this would be to divide the excess energy over the lev- 
els using equipartion, so 1/3 translational, 1/3 rotational and 1/3 
vibrational excitation, with a Boltzmann distribution for the sub- 
levels. The formation of molecular hydrogen, H2, in the ISM 
has been a long-standing problem. At extremely low metallic- 
ities, when no dust is prese nt, H2 either forms through the H ~ 
route (H + H ^ H2 + e ': lGlovedl20()l iLaunav et al.lll991h . 
or through a three-body reaction (H -H H H- H H2 -H H). 
The three-body reaction dominates at very high densities ( 
n > 10'' cm-^: iGalli & Pallal Il998h . For normal ISM condi- 
tions, H2 is predominantly formed on dust grains, and much 
progress has been made over the past two decades, both exper- 
imentially as well as theoretically, in understanding this (much 



faster ) formation of H2 on the surfaces of gra ins (Pirronel lo et al 
19991: iKatzet all 1 1999^, 'Cazaux & Tielens' '2002; Zec ho et al 



2002; 



2006; 



Cazaux & T ielens 2004; Shaet al. 2005; Horneks er et al 
Morisset et a l. 2004; Cup pen et al.l i2010b, and references 



therein). In general, the understanding is that CO is formed in the 
gas-phase, but the origin of gas-phase water is not so clear-cut. It 
is possible to form water in the gas-phase through either a chain 
of neutral-neutral reactions containing a number of temperature 
barriers, or through ion-molecule reactions when the gas is mod- 
erately ionized (Xe ~ 10"^). The formation of water on interstel- 
lar dust can also be an efficient route in diffuse and dense clouds 
(ICuppen & Herbstll2007l : ICuppen et ani2010ah . Several routes to 
form water can be considered, involving successive hydrogena- 
tion of oxygen either with H2 or H. Because species can stay on 
the dust, they can repetitively collide with each other and react, 
even if there is an important barrier for the reaction to occur. 
In this case, dust grains provide a favorable place for improba- 
ble reactions (with high barriers) to occur (Cazaux et al. 201(2). 
Also, other routes to form water, in volving O2 (Mivauc hi et al.l 
I2OO8I; Iloppoloetal.ll2010h and O3 (iMokrane et al.. ,20091 ) have 
been investigated and can be important at high dust tempera- 
tures (Tdust > 20 - 30 K). The molecules present on the sur- 
face can be released in the gas phase through several mecha- 
nisms: desorption upon formation if the reaction is exothermic; 
photo-dissociation of surface species can lead to the desorption 
of products in the gas (i.e., photodissociation of H2O leads to 
OH and H in the gas phase. [Andersson et alJl20 06): evaporation 
of the surface when the dust temperature is high enough (this 
process depends on the binding energy of the species). If many 
water molecules are formed and stay on the dust, then ices can 
be made. The bonds that hold the water molecules together are 
very strong (creating water clusters). As a result, when one or 
more monolayers of ice are formed on the grain, shocks will be 
the on ly effective way to remove the ice from grains (e.g.. lDrain3 
Il995h . It is thus only possible to have efficient formation of gas- 
phase water with grains as catalyst, when the grains are (close 
to) bare. 

Objects such as Mrk 231 have very violent environments in 
the central regions around the super-massive black hole. This 
particular object has a derived 8 - 10 00 yum infrared luminosity 
L,R = 4.0 X 10'^ Lo dBoksenberg et al.i,1977.). A highly absorbe d 
power-law X-ray spectrum was observed bv lBraito et alJ (l2004t) . 
It is in these extreme environments, that we expect cloud-cloud 
collisions to cause C and J shocks, and strip the grains of their 
ice layers. The internal UV field created by X-rays from the ac- 
creting black hole does not al low to reform the ice. The ab un- 
dances of H2O (xh,o ~ 10"^. iGonzalez- Alfonso et a n i20Tol) m 
the warm molecular gas of Mrk 23 1 are high. Here we study the 
effect of OH and H2O formation on bare dust grains in environ- 
ments with strong X-ray radiation fields, in order to investigate 
whether this process would be able to significantly boost the to- 



tal production rate of these species i n these environments. Our 
work differs from iHollenbach et al.l (l2009h . since we consider 
an additional process to release molecular species into the gas 
phase, and also since our study only considers bare interstellar 
dust grains. Even though this paper focuses on an application 
to ULIRGs, it will have applications to other objects as well, 
such as X-ray irradiated disks around Young Stellar Objects. The 
manuscript is articulated in the following way: We first derive 
analytical expressions for the formation rates of OH and H2O on 
dust grains. Then we compare these rates to the neutral-neutral 
gas-phase formation processes. Finally, the effects of the grains 
surface chemistry on the XDR thermo-chemical structure are 
shown. Their implications to explain recent observations with 
the Herschel Space Observatory are also discussed. 

2. H2O and OH formation rates 

As mentioned earlier, gas-phase OH and H2O can occur through 
neutral-neutral reactions, or when the medium is moderately ion- 
ized, Xg > 10"^, through ion-molecule reactions. The neutral- 
neutral reactions OH -H H2 -> H2O H- H and O -H H2 -> OH + H 
are slow at low tem peratures, because they have activation barri- 
ers (Wagner & Graff .1987} . and only contribute at temperatures 
r > 250 - 300 K). The other route is H+ h- O ^ 0+ -h H, 0+ h- 
H2 ^ OH+ + H, OH+ H- H2 ^ H2O+ + H, H2O+ -H H2 ^ H3O+ 
+ H, followed by recombination HsO^ -H e" ^ H2O + H. 

In our current X-ray chemistry code, only the grain surface 
reaction is taken into account for the formation of H2. In order 
to study the effect of OH and H2O molecule formation on grain 
surfaces, similar reactions for these species will be added to the 
chemical network 



n(0)vondustO"eoH, 
n(0)vondustO'eH,o, 



(1) 
(2) 



dn{OH) 
dt 

dn{H20) 
dt 

where eoH and eHzO efficiencies for the formation of OH and 
H2O, respectively, and vq the velocity of oxygen atoms. The 
expression for the efficiencies will be determined in the next 
section, where we calculate what fraction of an accreted oxy- 
gen atom is converted in gas phase OH or H2O. Two different 
grain size distributions are considered, an MRN ("Math is et al.l 
1977) distribution and a Weingartner & Draine (2001) distribu- 
tion, with Hdustc/nCH) - 10"^' cm"- and 2.8-10"^' cm"^, respec- 
tively. The MRN distribution considers a minimum dust grain 
size of 50A, while Weingartner and Draine (W&D) consider 
much smaller grains, which are actually PAHs, and can be as 
small as 5A. The W&D grain size distribution, however, seems 
more appropriate to environments such as the central regions of 
Mrk231, where dust grains can break in smaller particles due to, 
e.g., shocks. A larger cross section yield more enrichment of the 
medium with water and OH. The efficiencies are derived in the 
next section. 

3. Grain surface chemistry 

In this study, we consider carbon grains (graphite and carbona- 
ceous). In order to describe the population of a certain species on 
a grain, we need to consider additional processes apart from the 
reactions that normally occur in the gas phase, such as neutral- 
neutral reactions. These are: accretion to, evaporation from, tun- 
neling and thermal hopping of atoms and molecules on the sur- 
face of dust grains. 
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Table 1. Binding energies (in K) for carbon grains 



Physisorption 

H 
O 

H2 

OH 

O2 

H2O 

O3 

HO2 

Ho 62 



550 Bergeron et al. (2008) 

1390 Bergeron et al. (2008) 

600 Akai & Saito (2003) 

1360 Cuppen & Herbst (2007) 

1440 Cuppen & Herbst (2007) 

2000 Cuppen & Herb_st (2007) 

2240 Lee et al. (2009) 

2160 Cuppen & Herbst (2007) 

2240 Cuppen & Herbst (2007) 



Chemisorption: 



H 



8500 Jelea et al. (2008 



Species arriving from the gas phase arrive at a random time 
and location on the dust surface. The arrival time depends on the 
accretion rate, which can be written as 



Ra 



nxVxcrSx 



(3) 



where nx and vx are the densities and velocities, respectively, of 
the species X, and cr is the cross-section of the dust particle. Sx 
is the sticking coefficient of the species with the dust, and we 
consider Sx - 1- The species X can go back to the gas phase 
through evaporation. This rate can be written as 



R 



evap 



Vi exp(-Ex/kT) 



(4) 



with Vi is the oscillator factor of atom in site i (which is of or- 



der 10 s for physisorbed atoms), and Ex the binding energy 
of species X (see Table [T] for the adopted values). The bind- 
ing energy Ex can be either weak (for atoms and molecules) or 
strong (only for atoms ). In this study we c onsider binding ener- 



tug y ' 

gies similar to those in lCazaux et alJ(l2010l) . There are two types 
of interactions between the atoms and the surface: a weak one, 
called physisorption (van der Waals interaction), and a strong 
one, called chemisorption (covalent bound). Atoms from the 
gas phase can easily access the physisorbed sites and become 
physisorbed atoms. These weakly bound atoms can travel on 
the surface at very low dust temperatures, meeting each other 
to form molecules. Moving of specie s on dust grains occcur s 
through tunneling or thermal hopping dCazaux & Tielensll2004l) . 
The mobility of a species X moving from physisorbed sites to 
physisorbed sites can be written as: 



a,,,,(X) = vP.,JX) 



(5) 



where Ppp(X) is the probability for the atom to move from a ph- 
ysisorbed to another physisorbed site. In this study, we consider 
a range of surface temperatures T > 15 K. At these temperatures, 
thermal hopping dominates and the probability for a species X 
to go from a physisorbed site to another physisorbed site can be 
written as: 



Ppp(X) = exp{-Epp(X)/kT), 



(6) 



where Epp(X) is the barrier between to physisorbed sites, that 
we assume to be 2/3 of the binding energy Ex (see Table [T]). 
Hydrogen atoms that are physisorbed can also enter chemisorbed 
sites through tunneling effects or thermal hopping. Because of 



the high barrier against chemisorption I Shal (l2002h . H atoms usu- 
ally tunnel through the barrier to populate chemisorbed site. The 
mobility for the H atoms can be written as: 



a. 



pc ~ ^Ppc-: 



(7) 



where Ppc(ii) is the probability for a physisorbed H atom to enter 
a chemisorbed site (Cazaux & Tielens 2004). Once two species 
meet in the same site, they can form a new species if the ac- 
tivation baiTier for the formation can be overcome. According 
to the exothermicity of the reaction, the newly formed species 
is directly released into the gas phase. T his probability , calle d 
/des in this work, has been discussed in ICazaux et alj (l2010l) . 
and is the most important route to release products to the gas- 
phase. Species present on the surface can also receive a UV 
photon and become photo-dissociated. In X-ray exposed en- 
vironments, an internal UV field is created by excitation of 
Lyman-Werner bands, and about 30 percent of the local ly ab- 
sorbed X-ray is converted into UV dMalonev et al.|[T996h . The 
photo-dissociation of H2O on the dust can lead to the release in 
the gas phase of the two products of the dissociation, OH and 
H. T his mechanism has a very small chance to occur (2 per- 
cent, lAndersson & van Dishoecd 12008 ). Although this process 
is much faster than thermal desorption (the grains are too cold 
to make thermal desorption efficient), it is not the main contrib- 
utor in delivering species to the gas-phase. On grains, where ice 
layers are present, the desorption upon formation is less promi- 
nent, because the species are more strongly bound to the sur- 
face. In those environments, photo-desorption is the dominant 
route to release species into the gas phase. More details about 
the adopted reaction rates can be found in Cazaux et al. (20ldl). 
As we are interested in the oxygen chemistry, we will consider 
H, H2, O, OH, O2, H2O, O3, HO2, and H2O2 as grain surface 
species in the chemical network. 

Two different numerical methods are performed to follow the 
formation of species on dust grain surfaces: (1) a Monte Carlo 
and (2) a rate equations method. MC simulations are performed 
at several fixed dust temperatures. These calculations are very 
time consuming, and are used in order to check for which con- 
ditions the chemistry occuring on the dust enter the stochastic 
regime. The rate equations method is used to simulate the chem- 
istry occuring onto dust for a wide range of dust and gas temper- 
atures, gas phase chemical composition, radiation fields. This 
method is orders of magnitude faster than MC simulations. 



3. 1 . Monte Carlo method 

We use step-by-step Monte Carlo simulations to follow the 
chemistry occurring on dust grains. The dust grains are divided 
into square lattices of 10x10 adsorption sites. Each site on the 
grid corresponds to a physisorbed site and a chemisorbed site, 
so that the grain can be seen as two superimposed grids. Species 
originating from the gas phase arrive at a random time and loca- 
tion on the dust surface. This arrival time depends on the rate at 
which gas species collide with the grain (Eq.[3]). When a species 
arrives at a point of the grid, it can become physisorbed, or, if 
its chemisorption states exists and its energy is high enough to 
cross the barrier against chemisorption, it becomes chemisorbed. 
In our model, however, because of the high barrier to access 
chemisorbed sites, H atoms mostly arrive from the gas phase 
in physisorbed sites. The species present on the surface may re- 
turn to the gas phase if they evaporate (Eq.lDl. The species that 
arrive at a location on the surface can move randomly across 
the surface by means of tunneling effects and thermal hopping 



4 



R. Meijerink et al.: Enhanced H2O formation through dust grain chemistry in X-ray exposed environments 



(Eq. |5]l. Species present on the surface can also receive a UV 
photon and become dissociated at the rates shown in table I A. l] 
Once a chemical species present on the dust moves, or adsorbs 
a photon, or meets another species, the next event that occurs 
to this species is determined as well as its time of occurrence. 
Therefore, the events that concern every species on the dust are 
ordered by time of occuiTence, and for each event that occurs, a 
next event for the concerned species is determined. In figure [1] 
we follow the fraction of oxygen that is released in the gas un- 
der the form of OH and H2O for environments 1, 2, 4 and 5 (see 
Table |2] for adopted parameters). The parameters of these envi- 
ronments represent typical ambient low and high densities and 
radiation field s in the ISM of active galax y centers, and are simi- 
lar to those in lMeiierink & SpaansI ( l2005h . The solid and dashed 
lines show our results for Tdust - 20 and 30 K respectively. For 
each environments, 50-60 % (30-40%) of the oxygen arriving 
on the dust goes back in the gas phase under the form of H2O 
(OH). At higher dust temperatures (~ 35 - 40 K), on the other 
hand, species coming on the dust evaporate very fast. Therefore, 
the efliciencies decrease dramatically since the grain surfaces are 
only sporadically covered by species. 

3.2. Rate equations method 

The rate equations method can be applied when there is a mini- 
mum coverage of 1 atom on the dust. This condition is satisfied 
if the grains are not too warm, that species can encounter each 
other before evaporating. By comparing our rate equations re- 
sults with Monte Carlo simulations, we will define a range of 
dust temperatures for which rate equations are valid. 

The rate equation method couples the different equations for 
the different species on the dust. The processes described in the 
previous subsection are dependent on the dust and gas temper- 
ature, UV field, and exo-thermicity of the reaction occurring on 
the dust. The equation for species X can thus be written as fol- 
lows: 



= «gas(^)vxA^. - n{XMY)app(X) + 

(1 - fdes)n(YMZ)appiY) - RevapM^) 

-Rphotxn{X) + Rpho,yn(Y) (8) 

where app{X) and app{Y) are the mobilities of the species X and 
Y (see Eq. |5]l. The different terms of this equation present (1) 
accretion of species X from the gas phase, (2) the formation of 
the new species that involve species X, (3) creation of the species 
X by the encountering of species Y and Z that stay on the surface, 
(4) evaporation of species X, (5) photo-dissociation of X on the 
dust, and (6) the creation of X by photo-dissociation of species 
Y. 

3.3. Efficiencies derived from rate equation metiiod 

The formation efficiency is defined as the fraction of oxygen 
/(O) that is released in the gas phase under the form of another 
species. We studied this in a few different environments, chosen 
such that they resemble the conditions as found in the XDR mod- 
els of Meijerink & Spaans (2005). The gas temperature is fixed 
at a constant temperature T - 100 K, and dust temperatures 
ranging from Tdust - 10 - 60 K. Four different radiation fields 
were considered at two fixed densities, n - 10^ and 10^^ cm"^, 
in both the atomic and molecular regime. The parameters used 
are summarized in Table |2] 



Table 2. Models parameters 



Env. 




jc(HI) 


X(H2) 


x(OI) 


Go 


1 




1 




3.4 X 10-* 


10",10',102,103 


2 




10-' 


0.45 






3 




10-2 


0.5 






4 


16" 


1 


10-3 




loMoMoMo* 


5 




10-' 


0.45 






6 




10-2 


0.5 







In Figs. |2] and [3] we show the obtained efficiencies for 
OH and H2O. The figures show that between dust temperatures 
T^dust ~ 15 - 35 K, the efficiencies for OH and H2O are roughly 
constant. In this regime most of the oxygen goes to H2O. For 
dust temperatures Tdust ^ 35 - 45 K, the total efficiency for 
forming molecules is rapidly decreasing, and most of the atomic 
oxygen will not react to another species, before it evaporates 
from the grain. Therefore, the dust surface chemistry can only 
be effective, when the dust temperatures are not too high. Recall 
that the way to form OH and H2O is through O -H H -> OH 
and OH -H H — > H2O. For slightly higher dust temperatures, 
the chance of an OH molecule to evaporate from the grain be- 
comes higher than that for meeting another H atom. Therefore 
at temperatures Tdust > 35 - 40 K, most of the oxygen atoms are 
incorporated into OH, and then eon is highest. 

Increasing the radiation field, has a similar effect as increas- 
ing the dust temperature, i.e., it results in a lower fHio/fOH ratio, 
which is caused by photo-dissociation of H2O on the surface. 
The strongest effect is seen in the environments 3 and 6, where 
the dominant hydrogen fraction is contained in H2. In this envi- 
ronment, the hydrogen accretion rates are lowest, and thus also 
have the lowest atomic surface density on grains. When most of 
the gas is in atomic form, the hydrogen accretion is such that 
at high dust temperatures most of the oxygen is converted into 
OH. However, when the gas become more and more molecu- 
lar, the atomic hydrogen accretion rate onto the dust decreases, 
and a larger fraction of atomic oxygen ends up in O2. In envi- 
ronment 3, this is even the most likely product at temperature 
T'dust > 35 K. 

In Figs. ID and |5] the efficiencies for the formation of O2, 
O3, OH2, and H2O2 are shown. The efficiencies for these species 
show a peak at dust temperatures Tdust = 35 -40 K. The max- 
imum efficiency for producing O2 ranges from ~ 10-^ in the 
atomic gas to ~ 10-' in molecular gas, and generally the efficien- 
cies for other species are smaller. At temperatures T > 40 K, O2 
efficiencies dominate over those for OH, but this is the regime 
where efficiencies for the formation of molecular species be- 
come very small. 

3.4. Validity of rate equation metiiod 

In order to estimate the uncertainties in the MC calculation, 
and to check for consistancy between the MC and rate equa- 
tion method, we did two additional MC simulations at dust tem- 
peratures of 20, 30, and 35 (low density) and 40 (high den- 
sity). After a time f of ~ 1000 - 2000 years (low density) and 
~ 5 - 10 years (high density), the efficiencies are converged and 
oscillate around a stable value. At Tdust - 20 and 30 K, the H2O 
efficiencies vary between 0.5 to 0.6, and those for OH between 
0.3 to 0.4, and thus very similar to the rate equation method. At 




0.01 



100 



1000 

time (years) 



0.01 



0.1 



10 



time (years) 



Fig. 1. Results of our Monte Carlo simulations on a 10x10 sites grain (30 A) for environments 1,2, 4 and 5. The solid lines represent 
our results for T^ust = 20 K, dashed lines for Tdust = 30 K and dotted lines Tdust = 35 K (left panel) 40 K (right panel). 



the higher dust temperatures (35 and 40 K), the values obtained 
in the MC method are lower The rate equations method is over- 
estimating the efficiencies, and we are entering the stochastic 
regime. 



4. Analytical expressions for OH and H2O formation 
efficiencies 



We only concentrate on dust temperature T^ust > 10 - 15 K. 
At very low dust temperatures, OH and H20 are still able to 
form at a smaller efficiency, but the chemistry is different then 
for a number of reasons: (1) ices can form on the surface (2) the 
chemistry is different, because the binding energies are different, 
and (3) the surface of the grain will be saturated with H2. 



The rate equations meffiod gives the populations of the different 
species on the dust surface, which allows us to derive simplified 
analytic expression for formation efficiencies of OH and H2O, 
that is easy to use in numerical codes. This means that the small- 
est set of reactions needed to reproduce the derived efficiencies 
from the rate equation calculation is determined. Although the 
analytic expression of the efficiencies is derived in Appendix lAl 
the dependencies are given here: 
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fdes 



„n(0)n(H)Q'H ^ /de.soH'"(03)n(H)o'H 



RacAO) 
/rfeiH,o"(OH)n(H)Q'H 
RacciO) ■ 



RacciO) 



(9) 

(10) 



The fractions fdeson^ /rfe.!oH" ™d /rfcsn,o are 0.36, 0.25, and 
0.15, respectively. These are the fractions of the molecules that 
are released in the gas phase upon formation due to the exother- 
micity of the reaction. RgcdO) is the accretion rate of oxygen, 
while n(0), n(03), and n(OH) are the dust surface densities of 
O, O3, and OH, respectively, in monolayers, e.g., 1 monolayer - 
100% coverage. These surface densities are listed in Appendix 
lAl an is the mobility of hydrogen atoms on the grain surface. 
The efficiency for the formation of water does not contain a term 
with accretion of OH from the gas phase, since the accretion of 
O followed by OH formation dominates over OH accretion. 

A very important outcome of the simulation is that the for- 
mation of OH on grains is dominated by O -(- H — > OH up to 
temperatures T ~ 30 - 40 K depending on the environment, and 
by O3 H- H ^ OH + O2 at higher temperatures. This implies that 
apart from expressions for «(0) and «(OH), also those for n(02) 
and «(03) are needed. For completeness, we also give the one 
for n(H20). The importance of the O3 H- H reactions is shown 
in Fig. |6] It is obvious from the figure that the OH formation 
efficiency would be underestimated without the consideration 
of the reaction with O 3. These results are not inconsistent with 
those obtained by lopp olo et al.l (l2008l) . where the formation of 
H2O is dominated by the reaction sequence: H H- O2 — > HO2, 
H + HO2 ^ H2O2 followed by H -H H2O2 ^ H2O + OH. We 
do consider the same reaction rate paths in our w ork, but in this 
study, molecules are formed on bare grains, while lloppo lo et afl 
(l2008h consider icy grains mantles. In our simulation a very 
significant part of the molecules are desorbed upon formation, 
while the cycle of photodissociation of H2O and reformation 
on the grains significantly enhances the amount of water re- 
leased in the gas phase. It does not require the complex route 
through H2O2. There is also another way to form H2O that in- 
volves O3. The reaction O + O2 — » O3, that releases only a very 
small part into the gas phase upon formation (only a few per- 
cent). So although the efficiency for the formation of O3 is low, 
there is a signific ant amount of O3 on the surface to react to 
OH. The work by iHollenbach etall (2009) is also very differ- 
ent from ours. They also treat a surface chemistry network, with 
H2, OH, O2, and some other different species, such as CH, CH2, 
and CH3. They do not include processes with O3 and H2O2, 
intermediate species in the processes leading to water forma- 
tion. However, the biggest difference between our models are 
the dominant processeses to deliver species back into the gas- 
phase which are thermal, photo, and cosmic -ray desorption in 
the 'Hollenbach et al. (2009) paper, as they are considering ices 
in that work, whereas this is desorption upon formation in our 
work. They consider a PDR environment, and therefore they al- 
low for the formation of strongly bound ice-layers. We assume 
that this does not occur because of the turbulent motions strip- 
ping the grains and the strong photodissociating radiation fields 
that are internally created by X-rays. 

The approximate analytical fits to the rate equation results 
are very good for environments 1, 2, 4 and 5. However, environ- 
ment 3 and 6, which are highly molecular, start to show devia- 
tions for dust temperatures Tdust > 30 K for the H2O. This is due 
to additional reactions on the grain surface, that we did include 
in the rate equation method, but not in the approximate fit. 
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Fig. 6. Efficiencies for environment 3, largely molecular, low 
density n = lO-' cm"-' gas. Go = lO'*: Result from rate equations 
simulation (solid line), with (dotted line) and without (crosses) 
inclusion of the reaction O3 + H — > OH -1- O2 in the analytical 
expression. 

The fits should not be applied at very low dust temperatures 
(T^dust) < 10 K, because the analytical fits do not match the MC 
and rate equation results. 

5. OH and H2O formation: gas versus dust 

In order to estimate the importance of OH and H2O formation 
on dust grains, a comparison is made to the gas phase neutral- 
neutral formation rates which are given by (UMIST database: 
iLe Teuff et aPIIOOOtlWoodall et akllfool : 



W = 3.14-10-'^(^pexp(-^) cm^s-' (11) 



^H,0,gas = 2.05 ■ 10- 



-12 / -* 

300 



1736 , 
exp(-— — ) cm s , (12) 



while the conservative dust phase reaction rates (using the cross 
sections from the MRN distribution, ndu.stc/n(H) = 10"^' cm"^), 
are given by: 



T 

^OH.dust = 3.9-10 eoH(7'dust) cm-' s 



T 

gas 



Vo.dust = 3.9- 10 eH2o(7'dust) cm^ s \ 



(13) 
(14) 



Depending on the environment (i.e., the ambient atomic and 
molecular abundances, gas and dust temperature, etc.), OH and 
H2O form preferably on the dust grains or in the gas phase: 

- Although the efficiency for formation of OH and H2O on 
dust is strongly dependent on the dust temperature Tdust, the 
variations with gas temperature Tg^s at a fixed dust tempera- 
ture Tdust are very modest. 
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Fig. 7. Comparison of dust and gas phase reaction rates atn - 10^'^ cm ^ for an environment with low (left, xh = 1, x^^ = 3 ■ 10 
and high molecular abundances (right, xu - 10"^, xh,=0.5) as a function of gas temperature for dust temperatures Tdust = 15, 30, 
45, 60, 75 and 90 K. The OH abundances are varied between xqh = 10 ^ 10 ^ and 10 at low a nd 10'^, 10 ^ and 10' ^ at high 
molecular abundances. This resembles the unshielded and shielded environment of model 4 in .Meiierink & SpaansI (l2005h . 



- The gas phase formation rates have large activation barriers 
of Tgas - 3150 and 1736 K for the formation of OH and 
H2O, respectively. These channels are essentially cut off be- 
low temperatures Tgas 5 200 K. 

- Formation on dust depends on abundances of O and H, and 
in the gas phase on H2, O, and OH. 

For illustration purposes, a comparison between the gas and dust 
phase formation rates are shown for two different regimes in 
Fig- 13 ^ mostly mostly atomic (left) and molecular (right) envi- 
ronment. We do not show the ion-molecule reaction rates, even 
though they are considered in our chemical network. These pro- 
cesses can be more efficient than the neutral reaction chain at 
low gas temperatures, which is a domain that we do not con- 
sider in this study. The abundu nces resemble the un s hielde d and 
shielded region of model 4 in iMeiierink & SpaansI ( l2005h with 
density n = 10^^ cm"^: The atomic environment has abundances 
jch = 1, xh, = 3 ■ lO-\ xo = 3 ■ 10-4, xoH = lO-\ lO-\ and lO"'' 
(left), and the molecular environment has abundances xh = 10"^, 
xh, = 0.5, xo = 2 ■ 10-4, and xqh = 10"^ 10"^ and lO^^ (right). 
We adopt a range of OH abundances, to show the effect on the 
gas phase rates. We should note here that the OSU astrochem- 
istry databas43 (maintained by Eric Herbst) claims a smaller rate 
at high temperatures T > 500 for the H2 + OH reaction, and is 
given by ^o.gas = 8.4 x 10"'^ exp(-1040/rgas). This would 
make the dust surface reaction even more important. 

The efficiencies eoH and eH,o are reasonably high, ~ 0.3 and 
~ 0.6, respectively, at dust temperatures Tdust < 40 K. As a re- 
sult, the H2O dust formation rate is always dominating the gas 
phase formation in the atomic environment, while the OH dust 
formation is more important at temperatures Tgas < 800 K. Even 
though OH and H2O formation efficiencies drop very fast for 
higher dust temperatures, molecule formation with dust as cat- 
alyst will always become dominant below a certain gas temper- 
ature. The break even point ranges from Tgas ~ 100 - 400 K 
depending on the dust temperature T^ust =45-90 K. We should 
keep in mind though that for Tdust > 35 - 40 K, the formation 
rates are overestimated in our model. 

This picture is different in a molecular environment. Even 
when the OH and H2O formation efficiencies are high (Tdust < 
40 K, gas phase formation is more efficient at gas temperatures 

' http://www.physics.ohio-state.edu/ eric/research.html 



Tg-is > 100 - 200 K. Where in the atomic environment, the OH 
and H2O gas phase formation is suppressed by the very low 
abundances of H2, the dust phase formation rate is suppressed 
due to the low accretion rate of atomic hydrogen, which is de- 
pleted by two orders of magnitude. In the molecular environ- 
ment, the formation of water in either dust or gas phase is inef- 
ficient through neutral-neutral reactions, and ion-molecule reac- 
tions will also contribute significantly or even dominate, through 
the chain 0+ -H H2 ^ OH+ + H, OH+ -H H2 ^ H2O+ + H, 
H2O+ -H H2 -> H3O+ + H, followed by recombination H^O^ + 
e ^ H2O + H. 



6. Application to X-ray environments 

The inner regions of active galaxies from, e.g., Mrk 231, which 
are highly exposed to X-rays, are likely to exhibit enhanced for- 
mation of H2O on dust grains: (i) Dust temperatures are expected 
to be moderate, Tdust ~ 20 - 40 K, as the heating of dust is less 
efficient than by UV; (ii) The dynamical timescales are not very 
long and shocks are expected to be present, and therefore dust 
particles are not expected to be or stay covered with ice; (iii) 
Shocks may even be able to break up grains into smaller par- 
ticles, therefore enlarging the effective surface of dust, where 
water is expected to form. 

In Fig. |8] the chemical abundances of H, H2, OH and H2O 
and temperatures of an XDR model is shown with density n = 
lO''''' cm"^. Fx - 160 erg cm"^ s"\ and Solar metallicity, Z = 
Zq. Three different cases for the formation of OH and H2O are 
considered: 

1. Gas phase only: The OH abundance is xqh ~ 3 x 10^^ 
at low column densities, rising to ~ 10"^ just before the 
H/H2 transition, leveling of to 10"^. The H2O abundance is 
xh,o ~ 10"'° at column densities < 10^^ cm"^. As OH, 
H2O becomes more and more abundant toward the H/H2 
transition, and reaches an abundance xh,o ~ 10"^ beyond 
A^H ~ 10^^^ cm"^ 

2. Gas and dust with an MRN distribution: Nothing changes for 
OH. As shown in Fig.|7l the OH gas phase formation is more 
important at gas temperatures T > 800 K, at all dust temper- 
atures. The abundance of H2O, however, is higher by a fac- 
tor 8, at low column densities A^h < 10^^ cm"^, and thus the 
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Fig. 8. X-ray irradiated cloud with density n - 10^'^ cm"'' and radiation field Fx - 160 erg cm"^ s {left) The abundances are 
showns for H, H2, OH, and H2O in three different cases: formation of OH and H2O (i) in gas phase only (ii) in gas phase and on 
dust with an MRN size distribution, and (iii) in gas phase and on dust with a Weingartner & Draine size distribution, {right) Gas and 
dust temperatures for the same models. Note that the temperature profiles for the gas only models and those where OH and H2O 
formation on dust with an MRN distribution are included are the same. 



integrated column of warm water {T > 200 K). At low tem- 
peratures (i.e., at high column densities, A^h <: 10^^ cm"^), 
nothing changes for OH and H2O, as the formation route is 
dominated by the ion-molecule reactions. 
3. Gas and dust with a Weingartner & Draine distribution: The 
H2 formation rate is now higher as well due to the higher 
surface area of the dust grains. As a result a small increase 
is seen in the H2 abundance at low column densities. At the 
same time, the H/H2 transition occurs at lower column den- 
sities, A^H ~ 2 X lO^-' cm"^, compared to A^h ~ 3 x lO^-* cm"^ 
compared to the gas phase model. As a result, the OH peaks 
also at a lower column density, and reaches a slightly higher 
peak abundance. The OH enhancement is due to a higher 
gas phase formation rate, as a result of the higher H2 abun- 
dance, not because of the higher formation rate on dust. The 
water abundance in the high temperature range is enhanced 
by more than an order of magnitude. As the H/H2 transi- 
tion occurs at smaller column densities, the rise of the water 
abundance also occurs at these smaller column densities. 

To summarize, we find that the warm H2O colunmn density 
is enhanced by apprimately an order of magnitude over a col- 
umn A^H ~ 5 ■ 10^^ cm"^, when including H2O formation on 
dust. At larger column densities, where temperatures decrease 
and the gas becomes molecular, the formation on dust grains 
does not impact the abundance of H2O. The OH abundance is 
only indirectly affected by the higher H2 form ation rate on dust 
when using the IWeingartner & Draind ( 1200 lb distribution with 
the higher dust cross section. 

7. Discussion & Conclusions 

The role and importance of OH and H2O formation on dust is 
investigated. We use both Monte Carlo as well as rate equa- 
tion simulations to determine the efficiency at which oxygen is 
converted into OH and H2O on dust grains. The rate equation 
method is very fast, and is used to calculate a grid of different 
environment, low and high density (n = 10^ and 10^'^ cm"^), 
different molecular fractions, radiation fields, and temperatures 
(see Figs. 12] and [3]). The Monte Carlo simulations are used as a 
consistency check, and as these simulations are fairly slow, effi- 
ciencies are only calculated for a limited number of dust temper- 



atures. At dust temperatures Tdust < 35 - 40 K, the Monte Carlo 
and rate equations method yield the same results, within the un- 
certainties. At higher dust temperatures, the efficiencies obtained 
by the two methods start to deviate: the rate equation method 
systematically obtains larger efficiences than the MC method. 
This is where the process of OH and H2O formation on dust en- 
ters the stochastic regime. We fitted analytical expression to the 
derived results from the rate equations, such that they are easy to 
implement in any chemistry code. We find that: 

- OH formation on grains and release to gas phase is dom- 
inated by O H- H OH (rdust < 35 - 40 K, and by 
O3 -H H — > OH + O2 (for higher dust temperatures). The 
high dust temperatures results in decreasing hydrogen sur- 
face densities, because of the increasing evaporation rates. 

- H2O formation on grains and release to gas phase is domi- 
nated by OH -I- H ^ H2O. 

- The fraction of OH and H2O that is released in the gas 
upon formation is of the order of 36% and 9% respectively. 
However, in strong UV radiation fields, water on the dust 
can be photo-dissociated, and then reform to water, which 
then gives another contribution to the gas phase. This loop of 
forming-dissociating-forming.. increases the effective frac- 
tion of H2O released in the gas by a factor 5. 

- A gas and dust phase formation rates comparison (see Fig. 
|7]) shows that formation on dust grains is important for dust 
temperatures lower than 40 K. 

We implemented the analytical expressions for the forma- 
tion rates in the Meijerink & Spaans (2005) XDR code, and cal- 
culated a model with parameters representative for the inner re- 
gions of active galaxies. We find that the additional processes 
enhances the integrated column density of warm water. This in- 
creased column helps to explain unusually strong water lines ob- 
served for Mrk 231, that has an accreting supermassive black 
hole with X-ray luminosity Lx ~ 10'^* erg/s. This ULIRG shows 
for example a very high H2O 2o2 -In/ CO J - ^-1 ratio. These 
ratios are not observed in typical starforming environments that 
are mainly exposed to UV only (P DR) en vironments, such as 
the one observed in the Orion Bar (iHabar t et al. 2010) or M82 
dPanuzzo et al.ll20T3) . where this ratio is typically an order of 
magnitude or more smaller than in Mrk 231, and as such these 
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very bright H2O lines might be typical for regions that are ex- 
posed to high X-ray fluxes, close to the accreting black hole in 
active galaxies. 
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Appendix A: OH and H2O formation efficiencies 

The formation efficiency with which gas phase OH and H2O are formed can be written as a balance between the desorption fraction 
from the grain upon formation, and the accretion rate of oxygen atoms onto the grain. On grains, two formation processes play a 
major role for OH in different dust temperature regimes, while for H2O, there is only one major formation route: 

O + H ^ OH Tdust<40K (A.l) 

O3 + H ^ OH + O2 Tdust > 40K (A.2) 

OH + H ^ H2O (A.3) 

Therefore, the efficiencies should be written as: 
fdesoHn(0)n(n)aH /<i«oH' "(03)n(H)aH 

60H = 1 I A. 4 1 

RacciO) RacciO) 

/d«Hjo"(OH)«(H)Q;H 

= ^) ^^-'^ 

n(H), n(0), n(Oj,), and n(OH) are the surface density of H, O, O3, and OH on the dust. These hydrogen and oxygen atoms are 
accreted from the gas phase at a rate: 

RacciH) = "Hi^Hi monolayers (A. 6) 

Ns 

RacciO) = — — monolayers i \ (A.7) 
where nm and hqi are the densities of atomic hydrogen and oxygen in the gas phase, vm and vqi the mean velocities of these atoms: 



Vffl ~ 1.45 • 10^ cm s-i, (A.8) 

VQI ~ 3.6-104^^^cms-\ (A.9) 

and^Vs = 4- = S-IO^^ cm the number of sites per cm . H atoms arrive on the dust surface in physisorbed sites, and then continue 
to a chemisorbed site or evaporate from the surface. As a result, the density on the surface can be written as 

n(H) = — (A.IO) 

The mobility for hydrogen atoms to enter into a chemisorbed site is mostly due to tunneling, which we can then approximate by 

Opc = 10'^ exp(-2apc ■^ ^"'^2^'' ) = 6.7 X 10"^, where ape and Ea are the width and the height of the barrier between physisorbed and 

chemisorbed sites. The adopted values are ape = 1.5 A and Ea = 2870 K. The hydrogen surface density can then be calculated as 
follows 

1.45-10-i«nHiV¥ 

n(H) = ^ " (A. 11) 

6.7 X 10-3 + 1012 exp(^) 

1 dust ' 

The formation of OH and H2O can be determined analytically from the rate equations, when the most important processes are 
considered. The rate equations of O, OH, and H2O, as well as O2 and O3 are considered here, since the efficiency of OH depends on 
the abundance of O3 as well. The abundance change as a function of time for these species are given by by the following expressions: 



dnjO) 

dt 
dnjOR) 

dt 

t/»(H20) 
dt 

— ^ = (1 - fdesoJniOfao + (1 - /;„^^)n(OH)n(0)aoH - n(02)n(H)aH - Rphoto^niOj) - Revapo^niOj) (A.15) 
= (1 - fdeso,M02)n(0)ao - n(03)n(H)aH - n(03)n(OH)aoH - Rphoto^niO^) (A.16) 



RacciO) - n(0)n(H)aH + n(OH)/?pto,<,„ - 2niO)^ao - R^aponiO) (A. 12) 

(1 - /rf,,vo„)«(0)«(H)aH + Rpho,u,oniii20) - Rph„,„„n{OH) - «(OH)«(H)aH - i?..«poH«(OH) (A. 13) 

(1 - /d.^,o)«(OH)n(H)aH - /?pto,„ „n(H20) - n(H20)n(H)aH (A.14) 
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Table A.l. Adopted photodissociation rates in the approximation 



Reaction 




Rate(s"') 


OH + photon — > 


+ H 


3.9(-10) 


O2 + photon — > 


+ 


7.9(-10) 


H2O + photon — » 


OH + H 


8.0(-10) 


O3 + photon — > 


O2+O 


1.9(-09) 



When we consider the surface densities of the species at steady state, the surface densities of the species become 

n(0) = RacAO) ^^^^^ 

niH)an + Revapo - Rpho,o„ + V2/?acc(0)ao 

«(OH) = r J (A. 18) 

JdesH^o n(H)aH nTOarn 

(1 - /rf..„,o)n(OH)n(H)arH 

n(H20) = ' (A. 19) 

n(H)au + Kphotn^o 

(1 - fdeso,)n{^f»o + (1 - )«(OH)«(0)aoH 

= ^„ (A.20) 

(1 - /d«o,)"(02)«(0)ao 

«(03) = ^ L^. (A.21) 

n(H)a;H + n(OH)QroH + Rphoto3 

FilUng in the equations for the various included reactions, the expressions become as follows: 



n{0) 



n(OH) = 



n(H20) 



niOi) 



n(03) 



RacciO) 

10-exp(^)„(H) . lO-exp(^) - ^^^^ . 7/?..(O)10-exp(^) 



0.15+ 



(A.22) 



n(H) exp(^) + exp(^) - o.i5.3.9.ioia:S)/^.,)/«(H) + " ^O'^'^-^O) exp(^) 
0.64«(O) 

0,., 3.90.10--°Go , '""^'^PC^ 
n(H)10'2exp(^) n(H)10i2 exp(^) 

0.64«(O)«(H) 



0.15n(H) + 3.9 ■ IO-22G0 exp(:^) + exp(^) 

^ dust -* dust 

0.85n(H)n(OH)10i2 exp(r^) 
n(H)10i2exp(^) + 8 ■ IQ-'OGo 
0.85n(H)n(OH)exp(^) 



M(H)exp(^) + 8- IO-22G0 



(A.24) 



0.975n(OH)n(O)10i2exp(zM) + 0Mn{O2)\0^^ 

^ dust ^ dust 

n(H)10i2exp(^) + 7.90 • IQ-iOGo + lO'^expC^) 

J dust dust 

0.975n(OH)n(O) exp(^) + 0.64n(O2) exp(^) 

-* dust ^ dust 

n(H) exp(iP) + 7.90 • IQ-^^Gq + exp(^) 

-* dust ^ dust 

0.977n(O)n(O2)10i2 exp(^) 

^dust 

n(H)10i2exp(^) + 10i2exp(=P)n(OH) + 1.8 • IQ-iOGo 

J- dust dust 

0.977n(O)n(O2)exp(^) 



(A.25) 



n(H) exp(=P) + exp(^)n(OH) + 1.8 • IQ-^^Go 



(A.26) 



